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A WEB-page containing materials of comparing exper- 
imental data and UrQMD model calculations has been 
designed. The page provides its user with a variety of 
tasks solved with the help of the model, accuracy and/or 
quality of experimental data description, and so on. The 
page can be useful for new experimental data analysis, 
or new experimental research planning. 

The UrQMD model is cited in more than 272 publi- 
cations. Only 44 of them present original calculations. 
Their main results on the model are presented on the 
page. 
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Monte Carlo event generators play a very important role in high energy 
physics. One can mark the following areas of their application: 

1. Pragmatics or practical tasks - development of new or upgrade of 
old experimental setups to study some processes/interactions, design 
of detectors, Monte Carlo simulation of the detector responses and 
so on. Event generators applied at these should have fast operation 
speed, stability of work, and a rough reproduction of previous exper- 
imental results. As an example, let us mention applications of the 
UrQMD model [1, 2] for the development of detectors for research on 
nucleus- nucleus interactions (CBM collaboration [3]), and detectors 
for investigation of antiproton-proton annihilations (PANDA collabo- 
ration [4]) at future GSI accelerators. The RQMD and HIJING models 
have been used for analogous purposes for RHIC experiments. The 
well known GEANT package [5] is widely used to simulate various 
installations. 

2. Analysis of new experimental data and new investigations 
planning. They include comparison of new data with previous data 
and model predictions. As a rule, the new data do not agree with 
model predictions, so some questions arise in this regard: whether all 
special features of the setup have been taken into account; whether 
they are free of methodical errors; whether the theoretical models are 
used correctly; whether the model parameters set is right; whether the 
discrepancy between the experimental data and model predictions is 
of systematical character; whether the discrepancy has been observed 
in previous experiments; whether the discrepancy has been considered 
as an evidence of a new physical effect, and so on. They are solved 
differently, and often quite difficult using quite often event genera- 
tors. The generators should have a flexibility in parameters variation 
and physical scenario, as well as a sufficient physical meaning of the 
parameters. 

Another situation takes place at new research planning. A first 
question asked by experimentalists is related to the load of the setup 
by ordinary background processes. A second deals with the radiation 
condition of the experiment. A third question asks about the admix- 
ture of the background processes in the phenomenon under study, and 
how it can be damped, and so on. Clearly, experimentalists prefer to 
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use well approbated and well recommended models for their estima- 
tions. Here one can not avoid a study of model application experience. 
As a rule, there is no enough time for this. Thus, experimental col- 
laborations attract the authors of the models or use authors variants 
of model code to solve the questions. As an example, let us point out 
on estimations of secondary particles multiplicity in central gold-gold 
interactions obtained by CBM collaboration [6], and estimations of 
background processes intensities of the PANDA collaboration [7] . 

3. Scientific or cognition aims only - search for new effects or phe- 
nomena on the base of analysis of a discrepancy between experimental 
data and model predictions. One uses the fact that Monte Carlo mod- 
els are a synthesis of existing notions about process mechanics. Thus, 
the discovered discrepancy can be considered as an evidence of our 
insufficient understanding, or as an evidence of new effects. For ex- 
ample, the discrepancy between experimental data and intra-nuclear 
cascade model calculations growing with the collision energy rising led 
in its time to appearing a very important conception for high energy 
physics - " formation time of secondary particle" . 

The final aim of all the efforts is creation of a theory of processes that 
could predict effects with any predetermined exactness. As there are only 
few such theories, the aim is re-formulated - creation of a theory or a model 
predicting observable effects with specified exactness. Determination of the 
exactness is a special additional task. 

The philosophical aspect of the scientific research - "cognition of Good 
wisdom" , is out of the scope of our consideration. 

The aim of this work is to create a WEB-page containing materials for 
the second trend on the well-known Ultra-relativistic Quantum Molecular 
Dynamics model 2 ) application. This foresees the following: to compose a 
list of papers where the model is mentioned, to select papers with original 
results on the model, to classify materials and to create their graphical 
representation. 

According to the electronic database of scientific publications 
(http://www-spires.fnal.gov/spires/hep/search/), the milestone papers on 
the UrQMD model [1, 2] were mentioned in 272 publications at the be- 
ginning of year 2006. All of them were looked through, and some of them 
presenting materials on the model were selected. The last publications were 

2 A short description of the model see in Appendix 1 
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studied for original calculations with model usage. Only 44 publications 
were selected (the list of the papers is given in Appendix 2). Graphical 
materials from the papers were put on the WEB-page. 

All published calculations in comparison with experimental data were 
sorted according to the following sections. 

1. Production cross sections 

2. Particle multiplicities 

3. Particle multiplicity ratios 

4. Multiplicity distributions and correlations 

5. Rapidity distributions of pions, kaons, protons and others 

6. mT-distributions, temperature 

7. P^-distributions 

8. Flow 

9. Event-by- Event fluctuations 

10. Dileptons, J/Psi production 

11. Others 

Three subdivisions were introduced in each of the section - hadron- 
nucleon, hadron- nucleus, and nucleus- nucleus interactions. As a result, 
the main page looks as it is presented in Fig. 1. 

Each cell of the table contains references on the pages with correspond- 
ing materials. Each of the reference has a form XX_YY, where XX is an 
order number of the paper in the list of used papers, and YY is an order 
number of the figure in the paper. 

Section 1.1 (Production cross sections, hh-interactions) collects descrip- 
tions of cross sections of various reactions in hadron-nucleon collisions. 
Mainly, they have been presented in the first publications [1, 2]. Note that 
the set of the cross sections represents only a small part of cross sections 
collected in the well-known compilation [8] . One can think that an analysis 
of a larger set of experimental data allows one to define more precisely the 
boundary between the binary model of the hadron-hadron interactions im- 
plemented in the UrQMD model at low energies and the FRITIOF model 
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Figure 1: View of the main page. 



[9] used at high energies. This may improve the ratio between the cross 
sections of multi-meson and few-meson reactions in the UrQMD model. 

Section 2 (Particle multiplicities) gives calculations of multiplicities of 
7r-mesons, i^-mesons, protons, anti-protons and other particles. It is on 
the alert that the model underestimates the multiplicities of 7r-mesons, 
and overestimates the multiplicities of strange baryons, anti-protons in 
NN-interactions (see pages 1_3_11, 3_t89, 14_2). At the same time, the 
model overestimates meson multiplicities and underestimates mutiplicities 
of protons, anti-protons, and hyperons in nucleus-nucleus collisions (see 
Section 2.3). 

In Section 3 (Particle multiplicity ratios) relative multiplicities of K + /tt + , 
K /tt~ etc. are considered. They get actual due to the papers [10, 11, 12]. 
The problem with description of K + /tt + ratio in nucleus- nucleus interac- 
tions (see page 34_4) at energies larger than 10 GeV allows one to as- 
sume that at high energies an additional transverse pressure appears in 
nucleus- nucleus interactions [13] which is not implemented in the UrQMD 
1.3 model. 

In Section 4 (Multiplicity distributions and correlations) calculations 
of particle multiplicity distributions and correlation dependencies between 
particle multiplicities are collected. They can be of interest for experi- 
ments at sufficiently low energies with a good particle identification. Par- 
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tide multiplicity dependence on the impact parameter in nucleus-nucleus 
interactions, correlations between multiplicities of particles produced in 
different phase-space regions, correlations between energies of neutral and 
charged particles etc. can be put in the section. Corresponding calcula- 
tions are absent but the properties can be measured or have been measured 
in experiment. 

The main volume of results is presented in the Sec. 5 (Rapidity dis- 
tributions). Here i^ + -meson rapidity distribution in pp interactions at 
160 GeV is given which is in a bad agreement with experimental data (see 
page 23_2). At the same time, analogous theoretic calculations for nucleus- 
nucleus interactions agree with corresponding experimental data. Baryon 
rapidity distributions in ^-interactions as well as in nucleus-nucleus col- 
lisions at energies larger than 6 GeV/nucleon are badly described by the 
model, especially in the nuclear fragmentation regions (see 13_3, 23_1 and 
24_1). The worst situation takes place with the description of anti-proton 
production in P&P6-interactions at 160 GeV/nucleon. 

Section 6 (m^-distributions, temperature) deals with particle distribu- 
tions on transverse masses, %, and related estimations of particle source 
temperatures. According to the collected materials, meson spectra in pp 
and pA interactions are described quite well by the model. The problem ex- 
ists with A'-meson spectrum reproduction for nucleus-nucleus interactions 
at energies larger than 6 GeV/nucleon. One can think that one needs to 
take into account the creation of quark-gluon-plasma for its solution. 

In Section 7 (Py-distributions) not numerous calculations of properties 
of particle transverse momentum distributions are collected. The calcula- 
tions of photon transverse momentum distributions in P6P6-interactions 
are also presented here (see page 21_10). Since the existing version of the 
model (1.3) does not consider the so-called "hard" processes, one can not 
expect a good description of the distributions at Pt > 2 in the nearest 
future. 

Section 8 (Flow) presents calculations of characteristics of flow in nucleus- 
nucleus interactions (v\, v% p x , F). Taking into account that this trend 
gets quite actual after the RHIC result publication, the volume of corre- 
sponding model calculations does not look sufficient. In particular, there 
is no analysis of dependence on Pt or centrality of interactions. In order 
to recover this deficiency, we present calculations of ^(Pr) (pages 44_6 and 
44_7) without comparison with corresponding experimental data. 

Section 9 is devoted to Event-by- Event fluctuations. The characteristics 
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are expected to be actual in the nearest future. 

Section 10 collects characteristics of dileptons produced in hardon-nucleus 
and nucleus-nucleus interactions. 

Finally, various seldom calculations are presented in Sec. 11. 

The page is available at 
http:/ /hepweb.jinr.ru/urqmdl_3/validation/urqmdjnodel_validation.htm 

The authors are thankful to Prof. G.A. Ososkov and Prof. Kh. Abdel- 
Waged for a reading of the paper and important remarks. 

Appendix 1: Short description of the UrQMD model 

The Ultra-relativistic Quantum Molecular Dynamic model (UrQMD) [1, 2] 
is a microscopic model based on a phase space description of nuclear re- 
actions. It describes the phenomenology of hadronic interactions at low 
and intermediate energies (y/s < 5 GeV) in terms of interactions between 
known hadrons and their resonances. At higher energies, y/s > 5 GeV, 
the excitation of color strings and their subsequent fragmentation into 
hadrons are taking into account in the UrQMD model. The model was 
proposed mainly for a description of nucleus-nucleus interactions. Note, 
that up to now there is no unique theoretical description of the underlying 
hadron— hadron interactions, with their vastly different characteristics at 
different incident energies and in different kinematic intervals. Perturba- 
tive quantum chromodynamics (pQCD) can be applied to describe hard 
processes, i.e. processes with large four-momentum, Q 2 , transfer. But 
pQCD is formally inappropriate for the description of the soft interactions 
because of the absence of the large Q 2 — scale. Therefore, low— pt collisions 
are described in terms of phenomenological models. A vast variety of mod- 
els for hadronic- and nuclear collisions have been developed. The UrQMD 
model is the most appropriate one for energy range of the future PANDA 
and CBM experiments. 

The model is based on the covariant propagation of all hadrons consid- 
ered on the (quasi-)particle level on classical trajectories in combination 
with stochastic binary scatterings, color string formation and resonance de- 
cay. It represents a Monte Carlo solution of a large set of coupled partial 
integro-differential equations for the time evolution of the various phase 
space densities of particle species % = N, A, A, etc. The main ingredi- 
ents of the model are the cross sections of binary reactions, the two-body 
potentials and decay widths of resonances. 
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In the model, as in other Quantum Molecular Dynamics (QMD) models 
(see [15, 16, 17, 18]), each nucleon is represented by a coherent state of the 
form (h, c = 1) 

/ 2 \ 3 / 4 f 2 1 "I 

(f)i(x;qi,p h t) = I — J exp|--(f - q t {t)) 2 + ~m{t)x\ (1) 

which is characterized by 6 time-dependent parameters, q\ and pi, respec- 
tively. The parameter L, which is related to the extension of the wave 
packet in phase space, is fixed. The total n-body wave function is assumed 
to be the direct product of coherent states (1) 

$ = U(f>i{x,qi,Pi,t) (2) 

i 

The equations of motion of the many-body system is calculated by 
means of a generalized variational principle. The Hamiltonian H of the 
system contains a kinetic term and mutual interactions Vij (H = J2i T% + 
\HijVij). The time evolution of the parameters is obtained by the re- 
quirement that the action is stationary under the allowed variation of the 
wave function. This yields an Euler-Lagrange equation for each parameter. 

d(H) u d(H) 

P l = — and q l = — — . 3) 

dqi dpi 

#i = ^ + V A E«i>=V ft W (4) 
m j 

^i = -V fi E«i> = -VfiW (5) 

(Vij) = f d 3 xid 3 x 2 (f)*(f)*V(x 1 ,X2)(f)i(f)j (6) 

These are the time evolution equations which are solved numerically. The 
equations have the same structure as the classical Hamilton equations. 

The interaction is based on a non-relativistic density-dependent Skyrme- 
type equation of state with additional Yukawa- and Coulomb potentials. 
Momentum dependent potentials are not used - a Pauli-potential, however, 
may be included optionally. The Skyrme potential consists of a sum of 
two- and a three-body interaction terms. The two-body term, which has 
a linear density-dependence models the long range attractive component 
of the nucleon-nucleon interaction, whereas the three-body term with its 
quadratic density-dependence is responsible for the short range repulsive 
part of the interaction. The parameters of the components are connected 
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with the nuclear equation of state. Only the hard equation of state has 
been implemented into the current UrQMD model. 

A projectile or target nucleus is modeled according to the Fermi-gas 
ansatz. The wave-function of the nucleus is defined as the product wave- 
function of the single nucleon Gaussians. In configuration space the cen- 
troids of the Gaussians are randomly distributed within a sphere with 
radius R(A), 

m-r (\[A + (Ai-l)f r =(^-j\ (7) 

Po is the nuclear matter ground state density used in the UrQMD model. 

The phase-space density at the location of each nucleon is evaluated after 
its placement. If the phase-space density is too high (i.e. the respective 
area of the nucleus is already occupied by other nucleons), then the location 
of that nucleon is rejected and a new location is randomly chosen. 

The initial momenta of the nucleons are randomly chosen between and 
the local Thomas-Fermi-momentum: 

p™ ax = %c(?>TT 2 pf , (8) 

with p being the corresponding local proton- or neutron-density. 

A disadvantage of this type of initialization is that the initialized nuclei 
are not really in their ground-state with respect to the Hamiltonian used 
for their propagation. The parameters of the Hamiltonian were tuned to 
the equation of state of infinite nuclear matter and to properties of finite 
nuclei (such as their binding energy and their root mean square radius). 
If, however, the energy of the nucleons within the nucleus is minimized 
according to the Hamiltonian in a self-consistent fashion, then the nucleus 
would collapse to a single point in momentum space because the Pauli- 
principle has not been taken into account in the Hamiltonian. One can use 
a so-called Pauli-Potential [19] in the Hamiltonian. Its advantage is that 
the initialized nuclei remain absolutely stable whereas in the conventional 
initialization and propagation without the Pauli-Potential the nuclei start 
evaporating single nucleons after approximately 20 - 30 fm/c. A draw- 
back of the potential is that the kinetic momenta of the nucleons are not 
anymore equivalent to their canonic momenta, i.e. the nucleons carry the 
correct Fermi-momentum, but their velocity is zero. Furthermore, the 
Pauli-Potential leads to a wrong specific heat and changes the dynamics of 
fragment formation. 
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Impact parameter of a collision is sampled according to the quadratic 
measure (dW ~ bdb). At given impact parameter centers of projectile 
and target are placed along the collision axis in such a manner that a 
distance between surfaces of the projectile and the target is equal to 3 fm. 
Momenta of nucleons are transformed in the system where the projectile 
and target have equal velocities directed in different directions of the axis. 
After that the time propagation starts. During the calculation each particle 
is checked at the beginning of each time step whether it will collide within 
that time step. A collision between two hadrons will occur if d < ^a tot /7r, 
where d and a to t are the impact parameter of the hadrons and the total 
cross-section of the two hadrons, respectively. After each binary collision 
or decay the outgoing particles are checked for further collisions within the 
respective time step. 

In the UrQMD model the total cross-section a to t depends on the isospins 
of colliding particles, their flavor and the cm. energy. The total and 
elastic proton-proton and proton-neutron cross sections are well known 
[20]. Since their functional dependence on y/s shows at low energies a 
complicated shape, UrQMD uses a table-lookup for those cross sections. 
The neutron-neutron cross section is treated as equal to the proton-proton 
cross section (isospin-symmetry). In the high energy limit (y/s > 5 GeV) 
the CERN/HERA parameterization for the proton-proton cross section is 
used [20]. 

Baryon resonances are produced in two different ways, namely 

i) hard production: N+N^ AN, AA, N*N, etc. 

ii) soft production: 7r~+p^ A , K~+p^ A*... 

The formation of s-channel resonances is fitted to measured data. Partial 
cross-sections are used to calculate the relative weights for the different 
channels. 

There are six channels of the excitation of non-strange resonances in the 
UrQMD model, namely NN -> N A 1232 , N N* , N A* , A 1232 A 1232 , A 1232 N* , 
and A1232A*. The A1232 is explicitly listed, whereas higher excitations of 
the A resonance have been denoted as A*. For each of these 6 channels 
specific assumptions have been made with respect to the form of the ma- 
trix element, and the free parameters have been adjusted to the available 
experimental data. 

Meson-baryon cross-sections are dominated by the formation of s-channel 
resonances, i.e. the formation of a transient state of mass m = y/sTh, con- 
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taining the total cm. energy of the two incoming hadrons. On the quark 
level such a process implies that a quark from the baryon annihilates an 
antiquark from the incoming meson. Below 2.2 GeV cm. energy interme- 
diate resonance states get excited. At higher energies the quark-antiquark 
annihilation processes become less important. There, ^-channel excitations 
of the hadrons dominate, where the exchange of mesons and Pomeron ex- 
change determines the total cross-section of the MB interaction [21]. 

To describe the total meson-meson reaction cross-sections, the additive 
quark model and the principle of detailed balance, which assumes the re- 
versibility of the particle interactions are used. 

Resonance formation cross sections from the measured decay properties 
of the possible resonances up to cm. energies of 2.25 GeV/c 2 for baryon 
resonance and 1.7 GeV/c 2 in the case of MM and MB reactions have been 
calculated based on the principle. Above these energies collisions are mod- 
eled by the formation of s-channel string or, at higher energies (beginning 
at y/s = 3 GeV), by one/two ^-channel strings. In the strangeness channel 
elastic collisions are possible for those meson-baryon combinations which 
are not able to form a resonance, while the creation of ^-channel strings 
is always possible at sufficiently large energies. At high collision energies 
both cross section become equal due to quark counting rules. 

A parameterization proposed by Koch and Dover [22] is used in UrQMD 
model for baryon-antibaryon annihilation cross section. It is assumed 
that the antiproton-neutron annihilation cross section is identically to the 
antiproton-proton annihilation cross section. 

The total and elastic proton-antiproton cross sections are treated ac- 
cording to the CERN/HERA parameterization: 



with the laboratory-momentum p in GeV/c. The parameters are listed in 



For momenta pi a b < 5 GeV/c, UrQMD uses another parameterization 
to obtain a good fit to the data: 



a(p) = A + B p n + C ln 2 0) + D ln{p) 



(9) 



[2]- 




-2 



0.3 <p < 5 
p < 0.3 



(10) 




0.3 < p < 5 
p < 0.3 



(11) 
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For low lab-momenta the annihilation cross section is dominant. The 
sum of annihilation and elastic cross section, however, is smaller than the 
total cross section: 

Ad = (Jtot — CTel — CTann (12) 

This difference is called "diffractive" cross section in UrQMD, aam = Act, 
and is used to excite one (or both) of the collision partners to a resonance 
or to a string. In the string case the same excitation scheme as for proton- 
proton reactions is used. For high energies the "diffractive" cross section 
is the dominant contribution to the total antiproton-proton cross section. 

The final state of a baryon-antibaryon annihilation is generated via the 
formation of two meson-strings. The available cm. energy of the reaction 
is distributed in equal parts to the two strings which decay in the rest frame 
of the reaction. On the quark level this procedure implies the annihilation 
of a quark-antiquark pair and the reordering of the remaining constituent 
quarks into newly produced hadrons (additionally taking sea-quarks into 
account). This model for the baryon-antibaryon annihilation thus follows 
the topology of a rearrangement-graph. 

The hadron-hadron interactions at high energies are simulated in 3 
stages. According to the cross sections the type of interaction is defined: 
elastic, inelastic, antibaryon-baryon annihilation etc. In the case of inelas- 
tic collision with string excitation the kinematical characteristics of strings 
are determined. The strings between quark and diquark (antiquark) from 
the same hadron are produced. The strings have the continuous mass dis- 
tribution f(M) oc 1/M 2 with the masses M, limited by the total collision 
energy y/s: Mi + M2 < \fs. The rest of the yfs is equally distributed 
between the longitudinal momenta of two produced strings. 

The second stage of h-h interactions is connected with string fragmen- 
tation. The string break-up is treated iteratively: String — > hadron + 
smaller string. A quark-antiquark (or a diquark-antidiquark) pair is cre- 
ated and placed between leading constituent quark-antiquark (or diquark- 
quark) pair. Then a hadron is formed randomly on one of the end-points 
of the string. The quark content of the hadron determines its species and 
charge. In case of resonances the mass is determined according to a Breit- 
Wigner distribution. Finally, the energy-fraction of the string which is 
assigned to the newly created hadron is determined: After the hadron has 
been stochastically assigned a transverse momentum, the fraction of longi- 
tudinal momentum transferred from the string to the hadron is determined 
by the fragmentation function. The conservation laws are fulfilled. The 



12 



diquark is permitted to convert into mesons via the breaking of the diquark 
link. 

This iterative fragmentation process is repeated until the remaining en- 
ergy of the string gets too small for a further fragmentation. 

The fragmentation function f(x,rrit) represents the probability distri- 
bution for hadrons with the transverse mass m t to acquire the longitudinal 
momentum fraction x from the fragmenting string. One of the most com- 
mon fragmentation functions is the one used in the LUND model [23] . In 
UrQMD, different fragmentation functions are used for leading nucleons 
and newly produced particles, respectively: 

/ (x-B) 2 \ 

f(x) mic = exp I ^— p — 1 for leading nucleons (13) 

/(^)prod = (1 — x) 2 for produced particles (14) 

with A = 0.275 and B = 0.42. The fragmentation function /(x) pro d, used 
for newly produced particles, is the well-known Field-Feynman fragmenta- 
tion function [24, 25]. 

The fragmentation scheme determines formation time of created hadrons. 
Though there are various possibility (for details see Ref. [1, 2]). 

After the fragmentation, decay of the resonances proceeds according 
to the branching ratios compiled by the Particle Data Group [20]. The 
resonance decay products have isotropical distributions in the rest frame 
of the resonance. If a resonance is among the outgoing particles, its mass 
must first be determined according to a Breit-Wigner mass-distribution. If 
the resonance decays into N > 2 particles, then the corresponding N— body 
phase space is used to calculate their N momenta stochastically. 

The Pauli principle is applied to hadronic collisions or decays by blocking 
the final state if the outgoing phase space is occupied. 

The UrQMD collision term contains 55 different baryon species (in- 
cluding nucleon, delta and hyperon resonances with masses up to 2.25 
GeV/c 2 ) and 32 different meson species (including strange meson reso- 
nances), which are supplemented by their corresponding anti-particle and 
all isospin-projected states. The states can either be produced in string de- 
cays, s-channel collisions or resonance decays. For excitations with higher 
masses than 2 GeV/c 2 a string picture is used. Full baryon/antibaryon 
symmetry is included: The number of implemented baryons therefore de- 
fines the number of antibaryons in the model and the antibaryon-antibaryon 
interaction is defined via the baryon-baryon interaction cross sections. 
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Elementary cross sections are fitted to available proton-proton or pion- 
proton data. Isospin symmetry is used when possible in order to reduce 
the number of individual cross sections which have to be parameterized or 
tabulated. 

The UrQMD model reproduces nicely various properties of hadron- 
hadron and nucleus-nucleus interactions. 
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